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Suspended resonant optical cavities are basic building blocks for several experimental 
devices. An important issue is the control strategy required to bring them in the 
resonant or slightly detuned conhguration needed for their operation, the so-called 
locking procedure. This can be obtained with a feedback strategy, but the error 
signal needed is typically available only when the cavity is near the resonance with a 
precision of the order AL ~ where IF is the cavity hnesse and Xi is the laser’s 

wavelength. When the mirrors are freely swinging the locking can be attempted only 
in the short time windows when this condition is verihed. Typically this means that 
the procedure must be repeated several times, and that large forces must be applied. 
In this paper we describe a different strategy, which tries to take advantage by the 
fact that the dynamics of the mirrors is known at least in an approximate way. We 
argue that the locking procedure considered can be more efficient compared with 
the naive one, with a reduced needed maximal feedback. Finally we discuss possible 
generalizations and we point to future investigations. 
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I. INTRODUCTION 


Optical resonant cavities provide very effective ways to measure very small displacements, 
and are a key component of several experimental high precision devices, in particular interfer¬ 
ometric gravitational wave detectors such as Virgo AdvancecP and LIGO AdvancecP, which 
are currently near the start of new scientific runs with improved sensitivity, or KAGRA® 
which will follow in the near future. Other examples can be the macroscopic optomechanical 
devices proposed by several groupd^ for the production of squeezed stated using pondero- 
motive effects. 

In both cases the cavity mirrors must be able to move along the cavity axis, for gravita¬ 
tional wave detectors must couple efficiently to the tidal force of a gravitational wave and 
ponderomotive devices must respond to the radiation pressure force. This can be obtained 
by suspending the mirrors to a chain of pendula which allow for the longitudinal motion 
with the additional bonus of attenuating the seismic noise above some cut-off frequency f^. 

However seismic attenuation is not effective below fg, and the residual mirror motion 
could be too large to be compatible with a resonant cavity working in the linear regime: 
typically, for a fg of few Hertz the residual mean square displacement is of the order of 

10“®mpl. 

The cavity length must be stabilized with a control strategy operating below the cut off 
frequency but ineffective above it, where we need free motion. There are two problems that 
must be solved. The first is how to lock a cavity which is initially moving in a completely 
nonlinear regime, in a sense that will be precised later, with typical displacements for the 
mirrors of the order of few laser wavelengths. The second is how to maintain it in this 
linear (locked) regime once this is reached. In this paper we will deal only with the lock 
acquirement (locking) issue. 

In order to design a feedback control we need an error signal. This is provided by the 
phase shift experimented by a coherent light beam which enters the cavity and is reflected 
or transmitted. Let us consider a resonant cavity of length L built with two mirrors with 
reflection coefficient ri ^2 and transmission coefficient ti^ 2 - We are especially interested in the 
regime where both the reflectivities are near one. An important parameter is the cavity 
finesse iF, which is defined as the ratio between the full width at half-maximum bandwidth 
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of its resonances and its free spectral range and can be written as 

1 — rir2 1 — rir2 

As we will see in detail later it is possible to measure accurately the phase shift induced by 
the cavity on the reflected or transmitted light only near its resonance peak. For a suspended 
cavity which is freely moving the fraction of the time available for the measurement is thus 
given by and the typical duration of each passage in the measurement region will be 

of the order 


where v is the typical longitudinal speed of the mirrors, which depends mainly on seismic 
noise, and is the wavelength of the laser. If we want to stop the cavity during this time 
we need to apply a force such that 

mv _ mv"^ 

Fmax — ^ (3) 

r Xe 

which becomes large for an high hnesse cavity. Another point to consider is that the fre¬ 
quency bandwidth A/ of a feedback that must be present only during the measurement 
time r must satisfy 

A/ > — = ^ (4) 

47rr SttA^ 


In order to be able to design a simple and robust feedback the acceptable bandwidth must 
not extend too much, entering in a region where a large number of internal resonances of 
the mechanical system are present. Once again we can meet a problem when the hnesse (or 
the speed v) becomes large. 

Our proposal for a different approach to the problem starts from the consideration that 
by using measured data about the state of the system, only when they are available, we 
neglect a large amount of information which is known a priori. In a general framework 
we could assume that the state of the system is described by a hnite set of variables that 
we will denote collectively with x. We do not measure directly the state, but a different 
set of parameters y which are connected to the state in a probabilistic way, namely by a 
conditional probabilityP {y \ x) whose form is supposed to be known. We suppose also that 
we know the conditional probability of x at the time t given x' for t' < t, which we will 
denote P{x,t \ x',t'), and the initial a priori probability for the state P{x). 
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Our first aim is the estimation of the probability P {x, t \ yi, ti; ■ ■ ■ y^, tn) where t > tn, 
y^. is a previous measurement done at the time < tk+i and 0 < A; < n. Formally a simple 
recursive formula can be written for it using a prediction step which uses the knowledge 
about the dynamic of the system 


F (cc, f I fi; • ■ ■ ; y^, ^ J ^ ^ I ^ in \ yi,ti; ■ ■ ■ ;y^, tn) dxn (5) 

and an update step which integrates the information coming from a new measurement yn+i 

P (*n+l) in+1 I Vl) Fj ■ ■ ■ ) ?/n+l) ^n+l) — 


P {Vn+l 

^n+l) P tn+1 

- ■■ Wn^tn) 

I P {Vn+l 

^n+l) P {^n+litn+l 

Vli F) ■ ■ ■ ) Uni tn) dXn-\-l 


( 6 ) 


which is an application of the Bayes’ theorem. Our point is that the information contained 
in P [x, t I yi, ti, ■ ■ ■ ; yn, tn) is larger than the one inside P {x \ y), the last quantity being 
the only one that is used in the conventional approach. Our hope is that with this increased 
information it will be possible to improve the control strategy. 

The complicated multidimensional integrals involved in the previous equations can be 
computed explicitly only in very simple cases. For example, when the probability distri¬ 
butions involved are gaussian the previous steps lead to the Kalman hltei®®. Our locking 
strategy can be seen as inspired to some generalization of the Kalman hlter, as we will 
explain. 

In Section [IT] we will introduce a simplihed model for the system, an oscillator subject to 


white noise, which will be used to describe the locking problem in Section III Though the 
system basically can be described by a gaussian distribution in the state space, the problem 
of evaluating the locking performances reduces to the solution of a Fokker Planck equation 
with absorbing boundaries, which is not trivial and is studied numerically. 


We draw our conclusions in Section IV We propose some natural generalizations, and in 
particular a method that can be used to cope with poorly known parameters of the system, 
which is based on the update step of Eq. ([^. The technique can in principle be used to 
allow the identihcation of the suspension and the seismic noise model, but also of optical 
parameters. 
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II. THE MODEL 


The basic model we are interested to can be written in term of stochastic motion equations 
which we write in the Ito form 


dV + 2'yV dt + oJq {X dt — ag dWg) = ^Fe^t{t)dt 

dX = Vdt 


(7) 


Here X represents the variation of the length of the cavity with respect to a reference position 
Lq which is the equilibrium one when the laser is switched off. This is an harmonic oscillator 
of mass fi attached to a point which is subject to seismic motion in presence of a viscous 
damping proportional to 7 . The term Ff,xt is an external force which is introduced for future 
convenience. 

For the sake of simplicity the seismic motion is modeled as a white noise of spectral 
variance so dWg is a Wiener process, which is a quite crude approximation. A viscous 
modelization of the damping can also be inaccurate in some cases. 

More accurate models can be cumbersome in the time domain, though feasible. For 
example, in a suspended mirror dissipative effects are typically described by adding a small 
imaginary part to the stiffness constant k in the frequency domain (structural damping). 
These details would obscure the most relevant points we are interested to, so we ignore them 


here, however in the final section (Subsection IVA) we will discuss about a systematic way 
to introduce some improvement. 

A description of the cavity in the state space is easily obtained by grouping the dynamic 
variables as a; = Xei + UQ^Ve 2 . The stochastic motion equations become 


dx = F (x) dt + (3 dW 


( 8 ) 


where (setting f = /i Fg^t) 


Fix) = 


0 Wo 


0 


X 


= Ka; + /(f)e2 


-Wo -27 j \ fit) 

and (3 = Wodse 2 . The evolution of the probability distribution P [x, t \ x', t') is obtained by 
solving the Fokker-Planck equation induced by ([^ which reads 

dP {x, t I x', t') 


dt 


= -Vx ■ Jix,t) 


(9) 
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where 


J(a;, t) = F (a;) P {x,t \ x\ t') — -(3 {(3 ■ Va,) P {x,t \ x', t' 


( 10 ) 


can be interpreted as a probability current. The first contribution the the current in Equa¬ 


tion (10) gives a deterministic transport of the probability distribution, while the second 
one generates a diffusion effect proportional to the noise. Note that in our model the vector 
f3 has only the second component different from zero, and as a consequence the same is true 
also for the diffusion current. The reason is that the diffusion is generated by the random 
seismic noise force, which can randomize directly only the velocity. This randomization is 
next propagated to the position variable by the deterministic dynamic. 

The generating function Ws of the (connected) momenta of P over some region of interest 
S is defined by 

= [ P>-^P{x,t I x',t')(fx = 

Js 

and starting from the Fokker Planck equation (|^ it can be shown that it is governed by the 
equation 


(V,)e'^= + t(/3-i7)" (11) 

Here we neglected boundary terms which are not relevant if we are interested in the 
evolution of P over all the state space. By differentiating one and two times and setting 
77 = 0 we obtain the motion equation for the expectation value of the state vector 


m 


{x) 


{F) 


and the motion equation for its covariance array C = {x ^ x) — (x) ® (x) which is 


m 


c 


{x <Si F) + {F <Si x) - (x) (g) (F) - (F) ® (x) + f3<Si f3 


In our model F depends linearly on x, and we get 

k (j;) = K (x) + /(t)e 2 


( 12 ) 


which is exactly the equation of motion for the state without seismic noise and 


Ft 


c 


KC + CK^ + /3 0 /3 


(13) 
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which is unaffected by the external force /. 

The explicit solution of these equations can be written using the matrix U [t) = (see 


(14) 


Eq. (Al)). For the average value we get 

{x (t)) = U {t - t') {x {t')) + f U (t - r) 62 /(r)dr 

Jt' 

When the external force is absent (x) —)■ 0 on a time scale 7 “^, which is also the time scale 
in which the initial condition is forgotten. 

For the covariance array we obtain 


with 


C (t) = U (t - 1') C(T)U {t - t'f + Q{t- t') 


{t — t') = / U (t — r) {t — tY' dr 


(15) 


The initial value of C goes to zero exponentially on a timescale ( 27 ) ^ while the noise induces 


the contribution Q (see Equation (A2) in the Appendix for the explicit expression) which is 
initially zero and converges on the same time scale to 

ai 0 


lim Q = lim C = 

T^OO r^oo 


0 


with 


9 ~ 9 ~ 2 

2 ^ 

^00 
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which correspond to two uncorrelated components of the state vector with the same variance 
cr^. The parameter a'j introduced is the spectral variance of a force equivalent to the seismic 
displacement. 

The dynamic of our model preserve the gaussian character of a probability distribution, 
owing to its linearity, and a generic gaussian solution of the Fokker Planck equation can be 
writterfl^ 


P {x,t I x', t') = TV [x] X, C) (16) 

By setting (x (t')) = x' in Equation ([l4|), C {t') = 0 in Equation (15) and by substituting 


in Equation (16) we find in particular 


P{x,t I x', t') = Mx {t — t') x' + U {t — t) 62 f{T)dT, Q{t — t') ) (17) 

which correspond to the initial condition P {x, t' \ x', t') = 6 {x — x'). 
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Figure 1. The IPDHl signal x for a resonant cavity, as a function of the phase shift + Ldet)- 

In this particular case ri = 0.99, r 2 = 0.98 and Lq = 1™) which gives a free spectral range of 
fpSR = 150MHz. The modulation frequency is fm = (27r)“^a;m = 30MHz. The continuous line 
is the imaginary part xi = I™ X of tho IPDHl signal, the dotted line its real part XR = R'® X- The 
dotted hlled graph is the ratio r between the transmitted intensity and the input one. 

III. LOCKING 

The amplitude of the field reflected by the cavity can be written as When X 

changes slowly compared with the bandwidth of the cavity, which is a very good approxi¬ 
mation in the cases we are interested to, we can write 

^ r 2 e^fo — ri r 2 e^fo — ri 1 

1 — rir2e^fo 1 — rir2 1 -f- T (^1 — 

Here 0 = 27r (Lq -|- X) /Xi is the phase shift over the length of the cavity and the wave¬ 
length of the laser beam. We do not suppose that the cavity resonates at X = 0, so we 
write Lq = L^es + L^et where the resonant length L^es is an integer multiple of A ^/2 and 
— 1/4 < Ldet/Xi < 1/4. The reflectivities r*, the transmissivities ti and the losses of the 
mirrors are connected by the relation r? += 1 — with i = 1,2. 

The phase shift induced by the cavity can be measured with the Pound Drever Hall 
techniqud^^i^, which gives a complex signal x that can be written as 

x = XR + ixi = 'R (0) R- (0+)* - R- (0)* R- (0-) (18) 

and = c~^ {oj£ ± Um) (To + X) ~ 0±c“^a;mLo is the phase shift of the two [PDH] sidebands 
over the cavity length. An example of y is given in Figure [^for a somewhat arbitrary choice 
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Figure 2. A schematic representation of the locking procedure. Both the light transmitted by the 
cavity, measured by the photodiode PDT, and the phase shift Pound Drever signal measured by 
the photodiode PDR are used. The transmitted signal is used as a trigger, to discriminate between 
the inside ILRARI and the outside fLRARI regimes. Inside the ILRARl a standard feedback force FB 
is applied to the system, while outside the feed forward force FF is used. The feed forward force 
is provided by FT taking into account the value of the cavity velocity at the exit of ILRARl The 
low pass hlter allow for a smooth transition between feedback and feed forward. 


of the relevant parameters. It is evident that the connection between the cavity displacement 
and X is linear only in a small interval O aronnd the resonance, where it can be 

nsed as an error signal. 

Another nsefnl signal is given by the ratio between the transmitted intensity and the 
inpnt one, r = where and 


1 — rir2e2®‘^ 1 — rir2 1 + ^ (1 — * 


(19) 


In order for r to be different from zero we mnst allow for a (small) transmissivity ^ 2 - This is 
also plotted in Fignrej^ note that it can be nsed as a trigger, as it becomes different from 
zero in a signihcant way only in the linear region of the IPDHI signal. 

In Fignre we schematized the proposed improved locking strategy. The non dashed 
parts represent the nsnal locking scheme, the dashed one onr additions. The two parts 
are activated respectively inside and ontside the ILRARl in the hrst case we have a feed 
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backward scheme, in the second a feed forward one. 

We can assume that, up to a negligible measurement error, the state of the system will 
be known exiting the ILRARI The main problem for locking acquisition is the value of the 
velocity: if too large, the usual lock attempt will require a large feedback force. Our strategy 
can now be described as follows: when the cavity exit the ILRARl we will stop to measure its 
state. But, using our knowledge of the probability distribution during the out of resonance 
period, we apply a feed forward force / optimized in such a way to reduce the cavity speed, 
when it will reenter the ILRARI 

In particular, suppose the exiting velocity of the cavity to be v^xu > 0. The cavity 
can reenter the ILR ARI around the same resonance, with a hnal velocity n < 0, or around 
the hrst one on the right, with a hnal velocity n > 0. A key quantity will be the ILRARI 
reentering velocity distribution Prvd {v; /), which is a function of and a functional 

of /. Starting from Prvd {v; f) it will easy to evaluate the square several parameter of 
interest. 

In order to calculate Prvd f) need to solve the Fokker Planck equation with 

peculiar boundary conditions which represent absorbing barriers at the lLR A^ . In Figure 
we represented the strip of the state space between two ILR ARI of interest (double dashed 
lines) centered in xi = xr and xi = xr = xr + y- When the state of the cavity enters 
a ILRARI we want to stop to consider it, so we set P {x,t \ = 0 outside the strip. In 

our particular case there is not a diffusion current in the horizontal direction: this means 
that P needs not to be continuous across the vertical boundaries of the strip. However the 
probability current on the boundaries 

Ji {x, t) = P{x,t\ x', t') F (a;) ■ ^ 

= uoP {x,t I x',t') 

must be directed outward, and this gives the boundary conditions, valid for t > t'. 


and 


P {x,t I x', t') X 2 


xi=xl,X2>0 


0 


P{x,t I x',t')x 2 


X\=XB.,X2<0 


0 


which are represented by the black thick lines in Figure]^ 
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Figure 3. A sample of schematic trajectories which contribute to the ILRARI reentering velocity 
distribution Pjiv d j with the same initial velocity. The continuous line contributes to the n > 0 part 
of Prvd, the dashed and dotted lines to the n < 0 one. Note that the noise allows the cavity to 
invert several time its speed before reentering in the ILRAR.I (dotted line). For illustrative purposes 
the trajectories are discretized, and have been simulated as a sequence of a vertical step due to the 
noise and a transport one. The tick lines are the absorbing boundaries, where P{x,t \ XqPq) = 0. 


After finding the solntion P* which satisfy P* {x, t' \ x', t') = 6 (x — x') with x' = + 

ujQ^Vexit^2 we can calculate the positive velocity part of Prvd as 


-PwD 

JQ 


dt 


Xl=XR,X2=iOQ V 


POO 

= v P* 
Jo 


!! " i.« 

OJq V \0Jq Vexit 


dt (20) 


and the negative velocity part as 


Prvd f) — / Ji 


dt 


Xl=XL,X2=i^Q V 


POO 

= V P* 
Jo 


XL 




XL 

^0 "^exit 


,0 


dt (21) 


As a consequence of the absorbing boundaries P* is not a gaussian distribution, and 
its normalization is not conserved: its integral over the strip at the time t > t' gives the 
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probability for the cavity of not being reentered the ILRARI at that time. 

Supposing that we do not care about the particular ILR ARI where the system will be 
locked, we can design our feed forward force in such a way to maximize square velocity 
reduction probability 


Pred iyexiti /) 


^Pexit 


Prvd {v, V, 


2 
exit 


it, f) dv 


'^exit 


In the best case we could obtain Pred f) > this would guarantees the 

possibility of “cooling” systematically the cavity until the locking becomes possible. 

The maximization must be done over a class of functions / which satisfy some set of 
requirements: / must not be too large and with not a too large frequency bandwidth. We 
will write / as 

^+oo 

f{t)= / 


where |/o(t)| < To and K{t) represent a linear, time invariant, causal low pass hlter. 

Owing to the non trivial boundary conditions, we do not attempt here to hnd an analytical 
solution of the Fokker Planck equation, neither to determine the optimal control. Instead 
we evaluate numerically the performances of a set of simple control strategies: 


Strategy 1: When Vexu is large, the only option is to try to slow down the cavity until it 
enters the next ILRARI on the right. A possible strategy is to apply the most negative 
constant force available /o = —Tq. The real force will be given by 


fit) = -Fogit) (22) 

where git) = Kit — t')dt' is the unit step response of the hlter described by Kit). 
There are no free parameters. 

Strategy 2: A second possibility, which should be effective in the intermediate region for 
Vexit, is to initially accelerate the cavity and then to decelerate it. This would make 
possible to bring the cavity to the ILRARI on the right in a short time (reducing the 
effects of the diffusion) with a small hnal velocity. In this case the applied force will 
be 

fit) = Foigit)-2git-Ti)\ (23) 

with a free parameter ti to adjust. 
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Strategy 3: With the third strategy we attempt to bring back the cavity to its starting 
ILRARI This is expected to be effective for low enough values of Vexit- We need a 
deceleration phase followed by an acceleration one, needed to stop the typical cavity 
which is moving back. This gives 


/(t) = -Fa [g{t) - 2g{t - t,)] 


(24) 


and also in this case there is a single adjustable parameter. Note that the option 
can be seen as the Ti —)■ cxd limit of this one. 


A. Numerical results 

To be definite we choose a set of parameters which are somewhat representative of the 
typical scenario one encounter in an interferometric gravitational wave detector, setting 
ojq = 27rrads“^, 7 ~ 10“^a;o and = 10“®m. Setting an upper frequency cut-off around 
fc — lOOHz for both seismic noise and control force bandwidth we have also dg = 10“®/ ^/J~c — 
10“^mHz“^/^ (supposing the root mean square of seismic displacement to be 10“®m). We 
suppose the maximal force we can apply to the mirror to be of the order of 10“^N, which 
for a cavity mas^of /i = 20 kg gives To = 8 x 10 “®ms“^. 

It will be useful to give some order of magnitude estimates in absence of external forces. 
With the given parameters the length of the free cavity is spread over aoo — 10“®m, and its 
typical velocity is — 6 x 10 “®ms“^, which gives a typical time needed to move from a 

ILRARI to the next one of r-p ~ 4 x 10“^s. 

Note that both uqTt and 'jtt are small, so we expect the details of the dynamics to be 
relevant only for for velocities much smaller than the typical one. In the same typical regime 
the relative spread of the initial velocity can be approximated as a^/ (caoCToo) — 2 ^ 777 : once 
again we expect diffusion effects connected to the noise to be relevant only for velocities 
much smaller than the typical one. 

This means that when Vexu has a typical value, or larger, we can estimate the result of 
the feed forward force looking only at the average value of the probability distribution. In 
this case by applying the first strategy we obtain always a reduction of the cavity velocity. 
If we have an infinitely large frequency bandwidth at our disposal then f{t) = —Fq and the 
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final velocity will be 


V ~ 


''exit 


FqX^UJq — V^xit 


0.5 X 10“® 

^exit 


In the real case we should take into account the effect of the low pass filter impulse response 


K{t) in Equation (22): the force will need a time 0{u~^) to rise to the largest value available, 
and the velocity reduction will be suppressed if UcTt is small. In short, we are sure to obtain 
a high Pred, the velocity reduction can be small, but obviously always larger than the one 
obtainable in the ILR.ARI as tt is larger compared with the time spent inside the ILRARI by 
a factor P. 

What happens when Vexu is smaller than u^a^o is less obvious. In principle diffusion 
effect in the phase space due to the seismic noise could hamper the possibility of reducing 
the cavity speed. 

We evaluated Prvd f) by repeatedly integrating numerically the motion equa¬ 

tions (|^ with the appropriate initial position and velocity, using a simple leap frog scheme. 
This gives an ensemble of trajectories in the phase space (some of them are represented in 
Figure [^. When they reach the absorbing boundaries the integration is stopped and the 
final velocity is stored until a sufficient statistic is obtained. 

In Figure]^ the results for the first strategy are reported. We plot the cumulative proba¬ 
bility distribution of the ratio {v/vexu)'^, for a given Vexu which is given as a fraction of the 
typical one. The thick line correspond to the distribution obtained when the control force 
is present. For comparison the results without control force is also showed (thin line). 

The value of Fq is keep constant, and we applied a third order Butterworth low pass 
hlter. By looking at the value of the cumulative distribution when {v/vexof' = 1 we can see 
that we are able to obtain a velocity reduction probability larger than 1/2 when p = 10“^. 

When p = 5 X 10“^ and p = 1 we can’t appreciate the result, because as anticipated the 
final velocities are almost unchanged compared with the initial one. When p = 10“^ our 
feedback seems to obtain a results which is the opposite of the desired one. This can be 
understood because in this case the cavity has a velocity which is very small compared with 
the typical one: the noise fluctuations accelerate it, and when we apply the control force 
there is a larger time available for this acceleration. 

The general conclusion is that the first strategy seems to work in a specific range for Vexu- 
We did not attempt a full optimization at this stage. 
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Figure 4. Results for the cumulative probability distribution of the hnal square velocity when the 
hrst control strategy is used, for selected values of VexU = and Fq = 8 x 10“®ms“^ (bold 

lines). For comparison the results with Fq = 0 are also plotted (dotted lines). The value of the 
square velocity on the horizontal axis is normalized to the initial one. 


In Figure [^similar results are reported for the third strategy. We hxed a value for the free 
parameter ti which should be the optimal one in absence of fluctuations for p = 3 x 10“^. We 
see that we obtain the desired objective for p = 10“^ and p = 5 x 10“^, which is in agreement 
with the expectations. We did not attempt a full investigation in this case neither, but we 
expect that with an appropriate tuning of ri, and with an appropriate use of the second 
strategy when needed, it should be possible to slow down the cavity whatever its initial 
velocity will be. 
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Figure 5. Results for the cumulative probability distribution of the final square velocity when 
the third control strategy (with ti = 0.04s) is used, for selected values of Vexit = P^ 0 <^oo and 
Fo = 8 X 10 ®ms ^ (bold lines). For comparison the results with Fq = 0 are also plotted (dotted 
lines). The value of the square velocity on the horizontal axis is normalized to the initial one. 

IV. CONCLUSIONS AND PERSPECTIVES 

We proposed a simple generalization of the common scheme used for the lock acquisition, 
giving some initial numerical evidence that the generalized scheme could get better perfor¬ 
mances. The discussed strategy is “blind”, in the sense that the additional control force 
applied to the cavity when it is outside the resonance region is a feed forward one, designed 
using only the last known value of the state variables and the information about the cavity 
dynamics. Several details have been neglected in this paper: we aimed only to discuss the 
basic principles, and a detailed experimental understanding is needed to discover potentially 
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weak points. 

There are however a couple of improvements that can be foreseen, and will be the object 
of further investigation. 

A. Improving the model 

Our model of the cavity dynamics is quite simple. In a real situation, the mirrors are 
suspended to a complex attenuation system needed to reduce external seismic noise at the 
desired level. This means that the simple oscillator considered in this paper should be 
substituted by a chain of coupled ones. And, as mentioned initially, the model for the 
dissipation is a rough one. 

In a similar way, we modeled the seismic noise in the frequency band of interest as white 
noise, while in a realistic scenario it will have non trivial spectral peculiarities, namely it 
will be a colored gaussian noise process. 

We do not expect these neglected details to have a big impact on the results. As a matter 
of fact, during its permanence between the small region between two[LRAR^ it will be quite 
a good approximation to neglect completely the dependence of the mechanical force from 
the position. 

The specific dissipation model can have a larger impact: in the viscous case dissipation 
effects will be larger at higher velocity compared with structural ones, so the estimation of 
the efficiency of the feed forward strategy can be different. And in principle strong spectral 
peculiarities of the seismic noise, introducing time correlations, can make some difference. 

All these are modelization issues, that can lead to the introduction of some unknown 
parameters. We stress that also in our simplified model there are some parameters which 
are totally unknown (such as Ldet) or known with some uncertainty ( 7 , Wq)- 

In our initial discussion we introduced the prediction step described by Equation (|^ and 
the update step described by Equation (|^. In designing the feed forward scheme we used 
the prediction step only, but the update one play an important role when we need to cope 
with some unknown or partially known parameters p. The basic idea is to redefine the state 
variables writing 
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where the dimension of x can be larger than two to accommodate a more refined mechanical 
model (a suspension chain, a realistic damping mechanism). The motion equation for x will 
depend parametrically by the variables p: now we can consider the extended probability 
distribution P {x,p, t | ■ ■ ■) and write a Fokker Planck equation for it. 


A first possibility is to impose a trivial dynamic for p (no evolution at all) and to start 
from a given prior which describes our ignorance of the p’s values. After each prediction step 


the evolved P {x,p,t | ■ ■ ■) will be compared with a new measurement during the update 
step, which will select the regions of the probability space which better represent the real 
system. 

If needed we can add a diffusive dynamics for the variables p, with the possibility of 
adapting to slow drifts of the system. We mention that some of the parameters g p 

can be used to describe colored seismic noise This can be done by writing 



where F(t) is some parametrized filter function and dW a Wiener process. Once again we 
can model our knowledge of seismic noise with some prior, and we can introduce a drift for 
pseism allow the model to adapt. 

In our model we completely neglect the effect of radiation pressure. This is not a too bad 
approximation for our purposes, because radiation pressure effects are depressed outside the 
ILR ARl When the radiation pressure is large we expect however residual effect at the ILR ARI 
boundary, which can have an impact on the feed forward design. 

There are no problems in principle in introducing the radiation pressure in our model. 
An important difference is that there will not be a gaussian solution for the evolution of the 
cavity, as the equation of motion will be nonlinear. This is not a great complication, because 
also in the simplified model studied in this paper the probability distribution of interest is 
not gaussian, owing to the absorbing boundaries. 

We mention that in principle it could be possible to compensate radiation pressure effects 
by looking at the transmitted signal r, which is just proportional to the laser intensity inside 
the cavity. 
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B. Improving the locking strategy 


A further improvement in the locking procedure performances could be obtained unblind- 
ing (at least partially) the control strategy outside the ILR ARI region. This would convert 
our feed forward procedure in a feedback one. 

A possibility worth to be studied is the one of using the Kalman signal y (and the trans¬ 
mission signal r) in the nonlinear region. Here the problem is not the nonlinear dependence 
of y from the cavity position, but its non univocity. During the update step this leads to 
the generation of a non gaussian probability distribution. This can be parametrized with 
good accuracy as a gaussian misture, and we are led to the concept of particle filterJ^^^, 
which can be seen as a way to parametrize a generic (not necessarily gaussian) probability 
distribution in the state space. This has the advantage of a simple implementation which 
does not requires tricky linearizations, but comes at the expense of a larger computational 
cost. 

We are currently investigating all these issues, and we plan to report on them in a future 
paper. 
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Appendix A: Explicit expression of evolution operators 

We list for reference the explicit expressions of the operators which appears in the evolu¬ 
tion equations for our system. The fundamental quantity is the exponential of the matrix 
K 


U (t) = 


cos ^ sin ilt ^ sin ilt 


sinCt 


cos Qt — ^ sin Qt 


(Al) 


where = cUq — 7 ^. The inhomogeneous term of the covariance matrix is given by 


[T] = 


(1 - e--) I - ^ 


^7 „-27r 


sin f 2 tM(t) 


(A 2 ) 
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where 


^ sin Qt + cos Qt ^ sin Qt 
M(t) = ^ ^ 


y sin ilt ^ sin fit — cos fit J 

and I is the identity matrix. 

The explicit expression of a multivariate gaussian distribution for the variable x, with 
mean x and covariance matrix C is given by 


A/'(a;;a;,C) = ^—^L=exp -^{x-xfc ^ 


Appendix B: List of acronyms used 


PD: Probability Distribution 

PDH: Pound Drever Hall 

LRAR: Linear Region Around Resonance 
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